Dicke model: entanglement as a finite size effect 
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We analyze the Dicke model at zero temperature by matrix diagonalization to determine the 
entanglement in the ground state. In the infinite system limit the mean field approximation predicts 
a quantum phase transition from a non-interacting state to a Bose-Einstein condensate at a threshold 
coupling. We show that in a finite system the spin part of the ground state is a bipartite entangled 
state, which can be tested by probing two parts of the spin system separately, but only in a narrow 
regime around the threshold coupling. Around the resonance, the size of this regime is inversely 
proportional to the number of spins and shrinks down to zero for infinite systems. This spin 
entanglement is a non-perturbative effect and is also missed by the mean-field approximation. 
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Coherent interaction between electromagnetic photon 
fields and matter attracted interest a long time ago [l[ 
with renewed attention gained in the last decade due to 
significant developments in the experimental techniques 
in various areas of physics. Achievement of Bose Ein- 
stein condensation of cold atomic gases in electromag- 
netic traps enabled the coherent coupling of hyperfine 
states of 10 5 atoms to a single photon mode of an optical 
resonator [3] • Advances in the semiconductor technology 
allowed to obtain optical microcavities where electron- 
hole excitations inside the semiconductor quantum well 
are strongly coupled to an eigenmode of the optical res- 
onator [3j • Strong coupling of a single mode of a trans- 
mission line resonator to a Cooper pair box [4|] and a 
quantized mode of an optical crystal cavity to several 
semiconductor quantum dots Q have been demonstrated 
as a possible way to a quantum computing device 

Theoretical understanding of all these systems is based 
on a model proposed by Dicke pj which describes N spins 
1/2 (identical two-level systems) with splitting energy 
2e coupled to a single mode of electromagnetic field to. 
It was shown that this model is exactly diagaonalisable 
[§]. At zero temperature it undergoes a quantum phase 
transition from a non interacting state with unpopulated 
bosonic mode to a condensed state with with a highly 
populated bosonic mode f§] if coupling between the bo- 
son and a single spin g is greater than a threshold value. 
In the thermodynamic limit a phase transition occurs in 
the region of strong coupling if temperature is less than 
a critical temperature which can be described by a Bo- 
golyubov Hamiltonian similarly to the pairing model of 
superconductivity [IcJ • Recently, a variational wave func- 
tion approach to the generalised Dicke model was used 
[IH to describe Bose Einstein condensation of exciton- 
polaritons in a semiconductor optical cavity. 

In this paper we analyse the Dicke model at zero 
temperature for a finite N using matrix diagonalisa- 
tion methods. We find that for the particular coupling 
strength, 
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gion. The lower bound of the inequality is onset of the 
quantum transition described by the mean field theory 
MJ. The upper bound is the condition to have only 
singly populated bosonic mode in the ground state. The 
approximated value is the result of 1/N expansion around 
the resonance. 

In the thermodynamic limit the ground state pro- 
jected onto the subspace of N spins is not entangled 
as it is a product state. For a weak coupling below 
the quantum transition the ground state is a product 
of all unexcited single spin states. For a strong cou- 
pling above the transition threshold the ground state 
is also a product state of all single spin states as 
a result of the mean-field approximation. In a finite 
system the mean-field approximation is not applicable 
in a small region above the transition threshold where 
the expectation value of the boson is of the order of 
1 and its fluctuations are also of the order of 1. The 
ground state in this region is a superposition of the un- 
excited state of all spins and a spin state with only one 
single spin-flip excitation, the N spin W-state, \W) — 

(IU|...} + im...} + IUT ...} + ---)/ViV. The W- 
state can also be interpreted as 'magnon state' at van- 
ishing wave vector fl3| . 

Furthermore, the W-state can be considered as a bi- 
partite entangled state in the following sense. Dividing 
all spins into two groups the W-state is a Bell state 
in the subspace restricted by only one spin-flip excitation 
above the unexcited states of each group of the spins, see 
Eq. (fT3| . In the course of a bipartite measurement if 
the first group is found in the excited state then the sec- 
ond group is projected onto the unexcited state and if 
the first group is found in the unexcited state then the 
second group is projected onto the excited state. 

We diagonalize the Dicke model for N spins 1/2 cou- 
pled to a single bosonic mode 



ff = W 6t6 + e£^+^£(S+6 + ^6t) 
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where the sum runs over N spin- 1/2 operators Sj that 



obey the commutation relations [Sf,Sj] = e a p y dijk 
and b (tf ) is standard bosonic annihilation (creation) op- 
erator. 



the ground state of the spin subsystem is a bipartite en- 
tangled state and it is not entangled outside of this re- 
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The Dicke model possesses the following conserved 
quantities. One is the number of excitations of the cou- 
pled spin-boson system, 

L = n + J z , (3) 

expressed in terms of the z-component of total spin oper- 
ator J a = J2j , a = x,y, z, and the occupation number 
operator n — b^b of the boson mode. Note that J a and 
n are not conserved separately. The eigenvalues of L are 
the so-called cooperation numbers c, given by the sum 
of expectation values of n and J z . A second conserved 
quantity is the total spin, J 2 = J 2 + (J+J- + J- J+) /2 
with eigenvalues j (j + 1). 

We represent the Dicke Hamiltonian j2| in a basis 
where J 2 and L are block-diagonal. Within a block of 
given c and j the remaining degrees of freedom can be 
labeled by the eigenvalues m of J z . In the representation 
\c,j, m) each block has a tridiagonal form. The diagonal 
matrix elements represent the energies of states contain- 
ing c — m bosons and m excited spins, 

(m\H\m) — lu (c — m) + em. (4) 

The first above- and below-diagonal matrix elements are 
transition amplitudes connecting all pairs of states which 
differ by just one flipped spin, 

(m\H\m + 1) = i= \J (c — m) (j (j + 1) — m (m + 1)). 
v N 

(5) 

The size of each block is limited by the fact that —N/2 < 
m < N/2 for j — N/2. However ,the upper bound on m is 
further constrained by the cooperation number c which 
has the lower bound —N/2 for a block with no bosons 
present in the state. 

The ground state of Eq. ([2]) is the lowest energy state 
of all blocks Eqs. Il4l")l) . For different values of the pa- 
rameters to, e, and g the ground state can have different 
c's and j's (see below). 

We consider now the case with the bosonic (lu) and spin 
excitation energies (e) close to each other. For the un- 
coupled case, g = 0, the ground state contains no bosons 
and no excited spins (i.e. all spins are, say, down). From 
Eqs. |(3]) and Q we see that the cooperation number of 
this state is then c = —N/2 with corresponding ground 
state energy E = —Ne/2. For finite but still small cou- 
pling, g -C 1, we can use a perturbative approach to 
remove the boson mode via a Schrieffer- Wolff transfor- 
mation to obtain an effective spin-Hamiltonian [lfjj ]. 

where the spin-boson mixing is eliminated up to the sec- 
ond order in g [15], introducing an effective XY coupling 
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Figure 1: Cooperation number c of the ground state obtained 
by matrix diagonalization of the Dicke model for iV = 3. The 
first solid line from the left indicates the quantum transition 
to the strongly correlated non-perturbative regime, and the 
second solid line is determined by Eq. ([9]). The dashed line 
corresponds to the resonance e — lo. 



between the spins within a subband of given boson oc- 
cupation number n (with n being conserved under the 
effective Hamiltonian Eq. ([6])). Such a perturbative ap- 
proach gives a qualitatively correct description of the 
energies and wave functions until g crosses a threshold 
value g c where a transition to a strongly correlated non- 
perturbative regime takes place (see below). The cooper- 
ation number c of the ground state for the Dicke Hamilto- 
nian Eq. |(2|) is plotted in FigQ]for the solution for N = 3. 
There are several regimes: for small g the ground state of 
the system is defined by a regime (black) with c = —N/2 
and (n) — 0. Then, with increasing coupling a quan- 
tum phase transition at g c takes place to a new regime 
with c > —N/2 where (n) > 0, and where subbands with 
sharp bosonic occupation numbers no longer exist. 

Increasing g further above the transition threshold a 
sequence of states with c = 1 — N/2, c = 2 — N/2, and 
so on becomes subsequently the ground state due to the 
interaction energy. The ground state with c greater but 
close to —N/2 can not already be approximated with help 
from the perturbation theory as coupling is too strong 
and can not yet be approximated using the mean field 
approach as the fluctuations of the order parameter are 
too large compared to its mean. Here we diagonalize 
the Hamiltonian Eq. using the matrices Eqs. (|4l5p . 
Sizes of the matrices in this regime are limited by the 
cooperation number c. In the presentation |c, j, m) the 
matrices are not larger then 2x2 and 3x3 for c = 1 — N / 2 
and c = 2 — N/2 respectively. The lowest eigenenergies 
of these matrices are 
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for c = 1 - N/2 and 



E 2 = - 



(N-2)e 
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for c — 2 — N/2. The last expression is a result of ex- 
pansion at the resonance in powers of ui — e. Both E\ 
and E2 belong to subblocks with j = N/2 in accordance 
with a theorem from [T3|. Comparing Eo and E\ we 
find g c = \/uJe. The same condition to have a non zero 
population of the bosonic mode in the ground state was 
established in 0|. 

Comparing E\ and Ei near the resonance ui — e we 
find that a state with c = 2 — N/2 becomes the ground 
state when g exceeds some value 172 given by 
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e + u+ (ui-e) /(2N-1) 



2 [2^1 - 1/2N- 1 



(9) 



Thus, g < (72 together with g > g c define the regime 
of the model parameters where matrix diagonalization is 
the only way to study the Dicke model. The upper bound 
of Eq. (1]) coincides with Eq. © for e « to and TV > 1. 
Increasing g further, we can determine the boundaries 
between the ground states with different values of c ( i.e. 
c = 3 — N/2, 4 — N/2, . . .) by numerical diagonalization, 
as shown in FigQ]for N = 3. A result similar to Eqs. ([7j 
lUini) was obtained in but the descreete jumps in the 
cooperation number of the ground state were interpreted 
as an infinite sequence of instabilities. 

In the strong coupling regime, g ^ g c , the mean field 
approach provides a good approximation to the exact 
ground state. Indeed, introducing the expectation value 
of the bosonic operator B = (b) and neglecting quantum 
fluctuations around it, the Dicke Hamiltonian Eq. ([2]) 
becomes 



H = u)\B[ 
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(10) 



Eigenstates of this Hamiltonian are product states of the 
N spins, and thus are manifestly not entangled. Diago- 
nalization of the 2 x 2 - matrices for each spin and subse- 
quent minimization of the sum of the lowest eigenenergies 
over \B\ gives the following self-consistency (mean-field) 
equation 



(11) 



^4\B\ 2 g 2 /N + e 2 



which describes a quantum phase transition at a thresh- 
old value of the coupling strength g c that we have already 
found from the matrix diagonalization. 
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Figure 2: Expectation value of n calculated by matrix diago- 
nalization of Eqs. 114151) - solid line and the mean field result 
- dashed line. Plots are for N = 3 spins and at resonance 
e = ui. The grey area is the regime defined by Eq.(l) where 
the cooperation number of the ground state is c = 1 — N/2. 



The mean field approximation is satisfactory even for 
a system of only a few spins. To see this, we com- 
pare for N = 3 spins in Fig|2] the expectation value 
of n taken with respect to the exact ground state of 
Eqs. (141 5 1) with \B\ 2 given in Eq. l[TTj) at the resonance 
lo = e. There we see that the largest deviation is in the 
region of intermediate coupling strength where the coop- 
eration number of the exact ground state is c = 1 — N/2, 
grey area in FigFJ For weak coupling g < g c the 
ground state coincides with the non-interacting one and 
(n) = \B\ =0. The quantum fluctuations of n compared 
to its mean value are already small for the coupling cor- 
responding to the exact ground state with c = 4 — N/2, 

^(n 2 ) — (n) 2 ^j I (n) 2 ~ 0.2. Thus, in the strong coupling 

regime the approximation of neglecting these fluctuations 
in Eq. (fTTJ|l is already good for g > gi irrespective of N. 

The ground state can be characterized in terms 
of entanglement between different parts of the spin- 
subsystem. In the weak coupling regime, the ground 
state, being a direct product of (unexcited) spin states, 
is not entangled. In the regime of strong coupling, 
where the approximate Hamiltonian Eq. (fl0| is valid, 
the ground state is also a product of the individual spin 
states which, thus, also has no entanglement of any pair 
of spins. The ground state in the intermediate region 
Eq. (Q} has to be found by matrix diagonalization and 
will be analyzed below. 

Diagonalization of the matrix Eq. (|4l5p for c = 1 — N/2 
and g — ^/uJe gives a ground state in the regime defined 
by Eq. Jl} . We change the mixed spin-boson representa- 
tion |c, I, m) to a separate representation of spins and the 
boson \n,l,m), where n is the bosonic occupation num- 
ber. Then, tracing out the bosonic degree of freedom n 
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we obtain the reduced density matrix of the spins only, 
which in the representation of J z eigenstates \m) is given 

by 
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(12) 



There is a finite probability to find either all spins in the 

completely polarized state or in the W-state \—N/2 + 1). 

The W-state is a Bell state for N = 2. For N > 2 

it is a bipartite entangled state. The set of spins can 

be divided into two equal groups consisting of N/2 spins 

each (assuming N even). In the basis |mi) \rri2), mi and 

1(2) 

m,2 are the eigenvalues of the operators J z belonging 
to the first (second) group, the state \—N/2 + 1) is 
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i.e. a measurement of one group projects another group 
onto the definite state. Therefore in the intermediate 
region Eq.(l) the ground state of the Hamiltonian Eq. 
is a bipartite entangled state. Note that N spin W-state 
does not belong to the class of N spin entangled states 
prepared by squeezing The squeezing parameter is 



not defined for the N = 2 W-state and is greater than 
one for N > 2 W-states. 

Let the spins have different splitting energies ej in 
Eq.{2j instead of the case ej = e which we have 
considered in the paper. A small variation (ej) — 
ej <C 6j, (ej) being the average over splitting ener- 
gies of all spins, can be treated using perturbation the- 
ory and will not affect our results, Eqs.(l,13) much. 
For instance, the N = 2 W-state in Eq.(13) becomes 
((1 - (e x - e 2 )/V2g) \][) + ||t)) /V2 at the resonance 
e = u>. The amplitudes of different components in the sin- 
glet will alter slightly while the W-state will not change 
qualitatively. For a more detailed analysis of the inho- 
mogeneous Dicke model see [191 ] - 

In conclusion, we analyzed the Dicke model for a finite- 
size system by matrix diagonalization. We found that 
the ground state is a bipartite entangled state only in a 
narrow regime of parameters next to the quantum phase 
transition at zero temperature. The ground state in this 
regime cannot be obtained from mean-field theory which 
approximates the ground state as a product state. Also 
perturbation theory is not valid in this regime. For an 
infinite system width of the corresponding parametric re- 
gion vanishes as 1/N if system is close the resonance. 
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